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We consider relativistic superfluids where the C/(l) baryon symmetry is spontaneously broken. 
We show that all terms in the expansion of the quantum effective action for the Goldstone fleld can 
can be found to leading order in derivatives, once the equation of state is given. This enables one to 
find the all scattering amplitudes between Goldstone bosons at low energies. We apply this general 
result to quark matter at asymptotically high densities and derive the low-energy effective action to 
leading order in Qs, and show that terms containing the fifth and higher powers of the Goldstone 
field appear only at the ai level with known coefficients. 



Introduction. — The behavior of matter at high density 
is being actively explored in connection with the physics 
of neutron stars. One of the most important issues is su- 
perfluidity, which is predicted for both neutron liquid 0] 
and quark matter at densities sufHciently high so that u, 
d and s quarks are present at approximately equal num- 
bers. The latter phase, called that color-flavor locked 
(CFL) phase and the closely related ones with kaon 
or pion condensation, are the subjects of numerous recent 
investigations |, |,|,|,§. 

At length scales relevant for astrophysics, the most im- 
portant modes in the superfluids (neutron or quark) are 
the massless Goldstone bosons and the superfluid vor- 
tices. Naturally, one would like to be have an effective 
theory which contains only these degrees of freedom. It 
is in principle possible, at least for quark matter at very 
high densities, to derive such a theory from the micro- 
scopic theory (i.e., QCD): the decay constants and ve- 
locities of the Goldstone bosons has been computed in 
Ref. |, |, and the same techniques should be ex- 
tendable to nonlinear terms in the Lagrangian as well. 
However, the calculations are technically complex and 
becomes impossible at strong coupling. 

In this paper, we show that from the equation of state 
alone one can construct a fully nonlinear quantum ef- 
fective action that completely describes the low-energy 
dynamics of the superfluid Goldstone fleld. The result is 
given by Eq. (pO|); this effective action allow one to find, 
in particular, the scattering amplitude of an arbitrary 
number of Goldstone bosons to leading order in energy. 

Quantum effective action. — The quantum effective ac- 
tion is an useful device to treat the problem of spon- 
taneous symmetry breaking ||Tl[| . Consider an arbitrary 
field theory in which the ground state breaks a C/(l) sym- 
metry spontaneously, and let <I>(x) be an order parame- 
ter. We do not require $ to be an elementary field; for ex- 
ample, in the CFL phase the superfluid order parameter 
is a colorless six-quark composite operator. The nature 
of $ will not play any role in the subsequent discussion. 

The partition function is defined as 



Z[J] = / V(j)^ ex:p[iS + i d^x J[x)^x) 



(1) 



where the path integral is taken over all elementary fields 
(pi in the theory. The partition function is the generating 



functional for all Green functions constructed from $. 
On the other hand. 



W[J] = ~i\nZ[J] 



(2) 



generates the connected Feynman graphs constructed 
from $. The quantum effective action r[(f>] is defined 
as the Legendre transform of W^[J], 
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J = 



(5$ 



(3) 



In perturbation theory r[<I>] is the sum of one-particle 
irreducible graphs, with external legs cut off. In other 
words, r[<I>] reproduces, at tree level, all Green functions 
constructed from 

The configuration where r[$] attains its minimum is 
the ground-state average of $ at vanishing source J — 
0. In the case of spontaneous symmetry breaking, the 
minimum is degenerate. At the minima r[$] is equal to 
iy[0], and is proportional to the ground state energy, 



minr[$] = W[0\ = -EqT, 
{*} 



(4) 



where T is the total time interval, and Eq is the ground 
state energy of the Hamiltonian corresponding to the ac- 
tion S. 

Derivation of the effective action for the Goldstone 
mode. — We shall consider QCD at finite baryon chem- 
ical potential /i, which is described by the Lagrangian 



C = Co 



M- 
2^7 q, 



(5) 



where Cq is the usual QCD Lagrangian at zero chemical 
potential, q is the quark field. Note that Cq is Lorentz 
invariant, and the last term in Eq. (|^) is the only Lorentz- 
breaking term in the Lagrangian. This fact will become 
useful in our further discussion. The Lagrangian (|^) pos- 
sesses a U(l) symmetry generated by the baryon charge. 



(6) 



which is not broken by the QCD vacuum, but may be 
broken by the ground state at finite /i, say, due to color 
superconductivity Q. The order parameter $ ^ qqqqqq 
has baryon number 2. To be general, let us assume the 



2 



baryon charge of the order parameter is M; so U{1)b is 
broken to Zm- 

The entire formahsm of quantum effective action does 
not assumes Lorentz invariance at any stage and can be 
appUed to the theory (||). Because the Lagrangian con- 
tains ^ as an external parameter, F is a function of as 
well as a functional of 

r = rifi, $] . (7) 

For the sake of constraining the form of r[/i, $], it is 
useful to generalize Eq. (||) to 

C = Co + ^A^ix)qrq, (8) 

and treat as an arbitrary background field, and put 
Aq = fi, Ai ^ only at the last stage. Now T = r[^^, <i>], 
which is a functional of both the background A^^ and $. 
Since the Lagrangian (^) possesses a gauge symmetry 

q _> qe^<-/3 ^ A^^A^ + d^a , (9) 

the quantum effective action also respects this symmetry, 

r[A^, $] =r[A,,+5^a, $e*^^"]. (10) 

Our final goal is to derive the quantum effective action 
for the Goldstone field, which is the phase of <I>, 

so one should try to integrate out the amplitude |$| and 
keep only ip. In the effective action formalism, the pro- 
cedure of "integrating out" simply amounts to minimiza- 
tion with respect to |<i>|, because T should be used strictly 
at tree level. Therefore we define 

T[A,,^]^minr[A,,me^'^'^]. (12) 

{|*|} 

The gauge symmetry satisfied by r[Ap, is 

T[Af,, ^]=r[A,, + d^a, ^ + a]. (13) 

In particular, if a is constant in space-time, Eq. ( p^ ) 
implies the invariance of the effective action with respect 
to shifting ip by any constant value. 

Let us now make an expansion of r[^^, tp] in power 
series of the field (p as well as of the spatial derivatives. 
Due to the shifting symmetry mentioned above, at any 
given, say, n-th, power of ip the number of spatial deriva- 
tive is at least n: there should be no ip without a spatial 
derivative staying in front. If one keeps, at each power of 
ip, only terms with the smallest number of spatial deriva- 
tives (symbolically, only d"'(p" terms but not d™p'"' with 
m > n), then the effective action density £cff depends 
only on df^tp, 

T[A^, p] « fd^xCsiA^, d^ip) . (14) 



This will be the only approximation made to derive the 
final result (po[). After making this approximation, we 
still can use F to compute the scattering amplitude be- 
tween any number of Goldstone bosons to leading power 
in external momentum. For example, we can compute 
the leading ©(p") behavior of n-point Green functions of 
Goldstone bosons, but not the corrections suppressed by 
extra powers of momenta. In this way our effective action 
F becomes additionally "low-energy" in the sense of, e.g., 
the chiral perturbation theory, with the difference that 
the expansion parameter in our case is pip (momentum 
times field) instead of p (momentum) , and the effective 
action is strictly tree-level. We shall argue below that 
Eq. ( p^ ) is valid below the gap energy. 

The U(l) gauge symmetry ( [T^ further restricts the 
low-energy effective action to 

F[A^, ip] = Jd^xC^ffiD^p), (15) 

where D^ip = df^ip—Af^. Ccs now depends on one variable 
instead of two. In principle, £cff can depend on the field 
strenth F^,y — d^j^A^ — d^A^^ but since at the end we 
will substitute A^^ = (/i, 0), we can restrict ourselves to 
constant fields from now on. 

To find the form of £cff[-D^(^], we perform a minimiza- 
tion of F over ip. Assuming that the minimum is achieved 
at constant fields ip independent of a;, the result becomes 
V4£off [— ^/i] iVi is the total four- volume). If we can find 
this function, then Ccs[D^j^p\ is obtained by a simple re- 
placement —A^ D^ip. 

Because T[lp\ was obtained from F[$J by minimizing 
with respect to the modulus j^j [Eq. ([l^)], by taking the 
minimization over the phase ip we have performed a full 
minimization over all configurations of $. According to 
Eq. (^, Ccfs must concide, up to a sign, with the ground 
state energy density egg at constant external field A^j^, 

C,ii{-A^)^-eg_,{A^). (16) 

Moreover, since is the only Lorentz breaking term 
in the Lagrangian (||), the result must depend only on 
Af^A^. Therefore, it is sufficient to compute Cgg for A^^ 
aligned along the time axis. 

Let us look at egg for A^^ — (fi, 0) from the thermody- 
namic point of view. It is the ground state energy, per 
unit volume, of the Hamiltonian 

H = Ho- iiQb , (17) 

where is the vacuum Hamiltonian of QCD and is 
the total baryon number. We can rewrite that as egg — 
{{Hq) — I^{Qb))/V, where the averaging is taken over the 
ground state. We recognize that egg = —P, where P 
is the thermodynamic pressure at chemical potential /x. 
The dependence of P on /Lt cannot be further restricted; 
we shall assume that the equation of state P = P(/x) is 
given and move on. We now have 

CM'A^) = P {{A^An'^^) , (18) 
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and hence, the full low-energy effective action is com- 
pletely determined from the equation of state, 



1/2 



(19) 



Putting, in this expression, = (/i, 0), we obtain finally 



r[/i, = / d^xp[y/{doip~tiy-{d,ipy) 



(20) 



Again, the functional dependence of P in Eq. (po| ) on 
its argument is the same as in the equation of state P — 
P(/L^), which is an external element in our formalism. One 
can expand Eq. (|2^) over powers of dof and {diip)^ and 
find all vertices in the effective actions. 

Quark matter at high densities. — To illustrate how 
Eq. (|2^) can be applied a concrete situation, let us con- 
sider three-fiavor quark matter at very high chemical po- 
tential, where the strong coupling as is small and the 
ground state is the CFL superfluid (|]. The pressure, to 
leading order in as, is the same as for a free quark gas. 



p. ^ ^ciV/ 4 



(21) 



In the CFL phase Nc — Nf = 3; we howerver shall keep 
Nc and Nf explicit. The effective Lagrangian is now 



NcNf 



(22) 



(23) 
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It is interesting that in this case the series over (90) ter- 
minates at the four-power terms. The first term in the 
expansion (|2^) is the constant equilibrium pressure; the 
second term is a total derivative and can be neglected if 
one is not interested in vortices where (p is multi-valued 
(in which case it leads the Magnus force). The ratio of 
the coefficients of the two quadratic terms in Eq. ( |2^ ) 
gives the velocity of the Goldstone boson; as we see it is 
equal to l/-\/3, which is the hydrodynamic sound speed 
in an ultrarelativistic gas. 

Moreover, the relation determining the speed of the 
Goldstone boson Us, the pressure P and the energy den- 
sity e, 



dP 

97 



(24) 



is exact for any equation of state, not just for P ^ n*. To 
show that one expands Eq. 
and finds 



w to the quadratic order 



dP 



d^P 



n d^ 
fi dn 



(25) 



where n = dP/dfj, is the baryon density. By using the 
zero-temperature thermodynamic relations dP = ndn, 



de = ^dn, one arrives to Eq. (p4|). It is nontrivial that 
the speed of the Gondstone boson is exactly equal to 
that of the hydrodynamic sound, even as we are at zero 
temperature outside any hydrodynamic regime. 

The last two terms in Eqs. (03) describe the interac- 
tion between the Goldstone bosons. The knowledge of /i 
turns out to be sufficient to find the scattering amplitude 
between the Goldstone bosons. The Lagriangian ( [2^ ) 
should be useful for computing, e.g., the neutrino emis- 
sion and absorption rates JT2[ or the thermal conductiv- 
itity |l^ of the superfluid core of neutron stars. 

One can easily compute corrections to Eq. (^) by in- 
cluding corrections to the equation of state. Because at 
small coupling g the superconducting gap is exponen- 
tially suppressed as e~^/^ [Q, it never appears in the 
perturbative expansion of the pressure over as- One 
hence can use Ref . ||l^ , which provides P to the two-loop 
order. For example, with the first correction Eq. ( pT| ) be- 
comes 



NrN 
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1-^i^^^),^ (26) 

Nr TT 



The effective action is obtained by replacing /i in Eq. ( [2q ) 
by [(9^(/3— /i)^ — (9i(y9)^]^/^. One immediately sees that the 
running of as is essential to obtains terms of different 
structure than those appearing in Eq. (p3|); in particular, 
terms with five and more Goldstone fields appears only 
in the order with a coefficient proportional to the beta 
function coefficient (3a. 

Notice that if instead oi P ^ ji* the equation of state 
was P ^ fj,^ , then the superfluid mode would move with 
the speed of light and is non-interacting. Such a situation 
occurs, for example, in QCD at finite isospin chemical 
potential fij when ^ yU/ ^ '^t^/tt ||16|| . 

Hydrodynamic interpretation. — Although the effective 
action we have just derived is a fully quantum object, it 
is instructive to derive from it the classical equation of 
motion. We find 



dn{ no 



Mo 



= 0, 



(27) 



where the covariant derivatives are taken on the back- 
ground Af^ — (^4,0), and = {Dfi^ D^^Y^"^ and no = 
dP/dfi\^, 

=fj.a were introduced as notations. If we further- 
more denote 



(28) 



u --^ 
Mo 

so that M^M^i = I, then the equation of motion (|27 
obtains the form of the hydrodynamic conservation law 
(9^(nou^) = 0, where plays the role of the fluid ve- 
locity and no plays the role of the baryon density in the 
frame locally comoving with the fluid. 

Furthermore, let us construct the stress-energy tensor 
from the action (20). If we use the standard formula, 



(29) 
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then the conservation law d^T^^^ = is satisfied, but T^^ 
is not symmetric. One can correct the problem by adding 
to T^^'^ an extra term, 

T"" = f"" + A^tiqu^" , (30) 

which does not spoil the conservation law (since A'^ is 
constant and d^{nQU^) = 0), but renders T^'^ symmetric. 
Substituting Eq. (|^ into Eqs. (|2|) and (|3^), one finds 
that the new stress-energy tensor formally has the same 
form as that of a perfect fluid, 

j,^^ ^ a^Df'ipD"^ - 5^''£off = (e + P)u^'u•' - g"'' P . 
Mo 

(31) 

In order to perform the last transformation in Eq. ( |3l| ) 
one has to recall that CcS — P and use the thermo- 
dynamic relation noMo = e + -P, where both e and P 
are computed at the chemical potential /io- However, in 
contrast to the usual hydrodynamics where /io and 
are truly independent variables (except for the condition 
u^u^ = 1), in our case there is only one independent 
dynamical field, which can be chosen to be ip. 

To understand more fully the meaning of from the 
hydrodynamic point of view, let us go to the norelativis- 
tic limit (e.g., consider neutron fluid). This limit corre- 
sponds to baryon chemical potentials /j, is close to the 
nucleon mass to at. Assuming <C toat, then 

= —d,^ , (32) 

TOJV 

which is the well-known relation between the superfluid 
velocity and phase gradient in nonrelativistic superfluids 
[the mass staying in Eq. (|3^) is to^t instead of that for 
the Cooper pair 2m m because the factor of 2 has been 
absorbed in the definition of tp, Eq. (pi])]. 

The relativistic generalization of the superfluid hydro- 
dynamic equations has been discussed before p^ . 



Up to some field redefinitions, the equations derived 
in Refs. |l^ (at zero temperature) coincides with 
our equation of motion for ip. The approach taken in 
Refs. H, |l| is to start fr om the known nonrelativistic 
hydrodynamic equations and generalize them to the rel- 
ativistic case. We have arrived to the same equation fol- 
lowing a different path which emphasizes the role play by 
the symmetries (Lorentz and baryon number). We also 
achieved an understanding of nature of the action ( [20| ) 
as the low-energy quantum effective action for the Gold- 
stone mode, which can be used not only at the classical 
level but also for purely quantum calculations, for ex- 
ample, of the scattering amplitudes between Goldstone 
bosons. However, the formalism described here does not 
give the hydrodynamic equations at finite temperature. 

Conclusion. — It is surprising that so much information 
about the superfluid Goldstone bosons and their interac- 
tion can be infered from the equation of state alone. One 
should, however, excercise some care with the effective 
action (|2^). In this Lagragian terms with more deriva- 
tives than fields have been omitted, so its applicability is 
limited to energies below some scale. This scale must be 
of order of the gap energy A, which does not appear in 
the leading order effective action (po|). Indeed, explicit 
calculations |^ show that the dispersion curve for the 
Goldstone boson is no longer linear when its energy is 
comparable to the gap, which means that terms of the 
type p^if^ becomes important at that scale. It is natural 
to expect corrections to terms with higher powers of p 
to become important at the same scale. 
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by DOE Grant No. DOE-ER-41132 and by the Alfred 
P. Sloan Foundation. 
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